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O ! Abstract 
O 

j. ■ We present a novel approach for a systematic large-spin expansion of the t-J 

Hamiltonian which enables us to work without the constraint of no double 
^ ■ occupancy. In our scheme we can perform the large-spin limit ensuring that 

j-H ■ the low energy spin excitations are in exact correspondence with the physical 

o : ^ 

O ■ excitations of the s = ^ Hilbert space. As a consequence, we expect a smooth 



- 

X 



dependence of the physical quantities on the expansion parameter 1/s. As a 



■ first application of the method we study the case of a single hole in a Neel 

background. A systematic expansion in fluctuations about this stable solution 

indicates that by increasing t/J the quasiparticle weight strongly depends on 

the momentum carried by the hole. Results, obtained on small lattice sizes, 

are found in excellent agreement with exact diagonalization data. 
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The t-J model in two spatial dimensions is perhaps the most challenging "unsolved" 
problem in the theory of strongly correlated electrons, since it is now commonly ac- 
cepted to represent the low-energy Hamiltonian for the two dimensional copper-oxide high- 
temperature superconductors. Recent calculations based on the old fashioned, but reliable 
high-temperature expansion techniques, have indicated that spin-charge separation, obvi- 
ously present in this model in one dimension, may also characterize the elementary excita- 
tions in 2D, leading to a break-down of Fermi liquid theory and to a possible explanation of 
the anomalous properties of the high-temperature superconductors. [IH 

We consider Nh holes interacting by the t-J Hamiltonian 

HtJ = -tY. P{c{.C,,a + h.C.)P +JJ2iS^ S, - ^N, iV,), 

where {i,j) denotes a summation over the nearest neighbor (n.-n.) sites of the lattice, P 
is the projector onto the Hilbert space without doubly occupied sites. Si and Ni are the 
spin and number operators at site i, and cj^ and Cj o- are the usual creation and annihilation 
operators for electrons of spin cr. Henceforth we assume t > 0, J > 0. 

Due to the difficulty to deal with the projector of no double occupancy, several semi- 
classical approaches 0-|5[ leading to a mean-field description of Hfj have been proposed. 
In particular, in the large-spin approaches presented so far, the simplification of the t-J 
Hamiltonian is achieved by generalizing the spin-i polarization state of the electron to an 
arbitrary spin s state and by giving a fictitious spin (s — |) to the hole. The various methods 
differ for the definition of the enlarged Hamiltonian, a freedom left by the very fact that 
in a large-spin generalization one can only require that the physical Hilbert space, as well 
as the t-J model Htj, must be recovered for the value s = ^ of the expansion parameter. 
Unfortunately, all of them face the following fundamental difficulty: as soon as the extended 
Hilbert space is larger than the physical one, spurious low-energy elementary excitations 
emerge. Hence, it is not at all guaranteed that the low-lying excitations of the extended 
Hamiltonian correspond to some physical excitation of the original one, so that it becomes 
difficult or even impossible to derive reliable results by performing a systematic expansion 
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in fluctuations about the mean-fleld obtained by letting s — oo. 

Explicitly, in the Kane et al. model one obtains the Neel background as the mean-fleld 
solution for small J/t, but for s > | some of the allowed excitations change the spin of the 
hole, which is clearly unphysical. Other large-spin generalizations instead face the problem 
that the hole propagation becomes a nonperturbative process in the 1/s expansion. In the 
latter case, as was shown by two of us in Ref. one obtains phase separation even at small 
J/t, whereas it is now believed that the uniform ground-state is stable in the physical sector 

m- 

In this work we propose a new method to simplify at large spin the t-J model without 
facing the previous difficulties. In our approach we do not deal directly with if^j, but 
consider instead a natural extension of it as introduced by Sutherland , because the latter 
allows us to apply the spin-wave theory with a one-to-one correspondence with the physical 
excitations. 

In the t-J Hamiltonian the single site i can be occupied by 3 kinds of "objects": A hole 
(boson) |0)j, an electron of spin-up | t )i = 4,t|0)«' ^"^^ electron of spin-down = 
c|jO)j, (fermions) whereas, apart for an irrelevant energy shift, Htj can be thought of as 
the operator permuting pairs of n.-n. objects, with weight t for permutations of objects of 
opposite statistics, and weight J/2 (—J/2) for permutations of fermions (bosons). In order 
to work without the local constraint of no doubly occupancy, we consider the extended 
Hamiltonian H acting on objects of two fermion and two boson species by permutation of 
pairs of neighboring objects with the same weights as in the t-J model. Because the number 
of objects of a given species is conserved by construction, the reduction to the physical 
model can be obtained by projecting onto the invariant subspace where one boson species is 
absent. Hence, in our approach the projector operator P is washed out, as for the projection 
amounts just to flx a conserved quantity. 

To represent the extended Hamiltonian iJ in a way suited for our developments, we 
denote the fermion and boson objects at site i with the symbols |lcr) and |Ocr) {a =t)i); 
respectively, and use the representation |0 t)j = |1 T)j = l"^)? |0 i )i = fiQi-\'^)i 
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|1 i)i = Qi.-I^)) where // and Qi- are a spinless fermion creation operator and a spin-| 
lowering operator, respectively. The extended Hamiltonian reads 

H=Y.\-i + fjfhx^, + ^(1 - - n,)(x,, - 1) , 

{id) ^ ^ ^ 

where = 2(5i Qj + | and rij = ///j. A noticeable feature of the proposed model if, to 
be contrasted with the case of Htj in a slave-fermion representation [^Q, is the presence of 
the permutator operator x^j- both in the magnetic and kinetic part, as well as the bilinear 
dependence on the fermion operators. 

Henceforth we preserve the name of electron and hole for the two "particles" |1(t) and 
I Oct), respectively, and introduce the two commuting vector operators satisfying the algebra 
of the angular momentum 

5 = 5](l-n,)Q„ L = 5]n,4, (1) 

i i 

The operators S and L act nontrivially on \ \a) and |0(t), respectively, and accordingly S will 
be referred to as the physical spin and L as the pseudospin. The analogy with the properties 
of the spin and pseudospin operators is not only formal and will be discussed elsewhere |]TU 



In the following we shall also refer to Q = L ^ S as, the isospin vector. The operators 
(HI) commute with if, so that the quantum numbers spin S^-, total spin 5, and pseudospin 
L^, total pseudospin L associated to both "particles" are conserved. The physical Hilbert 
space of Htj corresponds to the sector where the pseudospin attains its maximum value 
= L = iVj2. 

Sutherland [0 has shown quite generally that the ground-state of the Hamiltonian {H) 
is at most degenerate with the physical one with maximum pseudospin L^. Unfortunately, 
this statement is rigorously valid only in one dimension, and for the case of a single hole in 
any dimension. The latter case is of course trivial, because for one hole the pseudospin is by 
definition equal to the maximum value L = ^. The proof presented in Ref. is not valid 
in 2D. In fact, following the reasoning one would obtain that for J = the ground-state 
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of Ht,j is the fully polarized Nagaoka state, whereas it is known that for large doping the 
singlet Gutzwiller projected Fermi gas has macroscopically lower energy ||11|| . 

However, the Sutherland's result is true for the special case t = and probably remains 
valid for physically acceptable J/t > in D>2. Hence, we shall leave it as a "conjecture". 
The importance of this conjecture is easily understood by noting that whenever it is satisfied, 
one can evaluate ground-state properties avoiding even the projection onto the Lz = Nh/2 
sector. 

The Hamiltonian H still represents a highly nontrivial problem and we now consider the 
large-spin approach allowing to simplify the model. Noting that Q is a irreducible spin-^ 
operator, we consider arbitrary higher-dimensional representations of the isospin vector and 
define the enlarged Hamiltonian Hg by substituting in the extended Hamiltonian H the 
permutator operator Xi.j with the rotationally invariant expression 

x^,^xS = ^Q^Q^ + l- (2) 

The overall factors and constants - irrelevant in the undoped case - are set by the require- 
ment that (Xij) is one or zero if the isospins of the particles at sites i,j are parallel or 
antiparallel, respectively. In our approach the Hilbert space is generalized by giving a ficti- 
tious spin-s both to the electron (i.e., |lcr)) and to the hole (i.e., |Ocr)) and in this respect it 
is quite different from the large-spin approaches proposed so far. 

Because at zero doping isospin and physical spin coincide (Z = 0), our approach leads to 
the conventional spin wave-expansion for the Heisenberg antiferromagnet: The ground-state 
is a singlet and the physical spin-wave excitations have S = 1, i.e., they are independent of 
the magnitude of the spin s of the extended Hilbert space. We see that in the large-spin 
limit the Hilbert space at each single site grows with s but the low-energy excitations remain 
in one-to-one correspondence with those of the physical s = | Hilbert space. We believe 
that this is the basic reason why the spin-wave expansion is so accurate for the undoped 
system and why 1/s is a smooth parameter and actually small i^Jl^-[l4 . 



At nonzero doping, for s > ^ the observable (|l]) are no more conserved, unless for t = 0. 
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However the Hamiltonian Hg always commutes with the isospin and it is therefore convenient 
to work in the subspace with Q and Qz fixed. If for the original t-J Hamiltonian the physical 
spin attains its minimum value - as it is expected for J not too small - the Hilbert space 
with minimum S and fixed 5*2 is exactly equivalent to the one with Q = [Nh/2] and say 
Qz = (for Nh odd one has S = 1/2 as the minimum value). Hence, if we are able to classify 
the elementary excitations for s ^ oo in the sector where L is frozen to its maximum value, 
the same classification would hold in the physical Hilbert space, provided that there is no 
phase transition function of 1/s. 

In order to show that our 1/s expansion is the natural extension of the spin-wave theory 
even in the doped case, we focus our attention to the simplest possible nontrivial doping, 
i.e., when there is only one hole {N^ = 1) in a lattice of M sites with periodic boundary 
conditions. 

For large s and J not too small the stable semiclassical solution corresponds to a Neel 
background where the single hole can propagate with given momentum either on the A or 
on the B sublattice. Fluctuations over this semiclassical solution are obtained in the usual 
way by introducing boson operators a| ~ ii i G A ii i E B) that create a spin 

fiuctuation over the Neel classical state \N). Thus a systematic expansion of the operator 
Xif in 1/s is possible and we get 



where ipij = a| + aj. By replacing the expression in the Hamiltonian H, both in the 
kinetic and the magnetic term, we then obtain an effective Hamiltonian for the single hole, 
which is characterized by a kinetic term coupling two boson and two fermion operators 
"2j(/j/l + h.c.)[(aj + aj){ai + a^) — 1]. Hence, it is remarkably different from the Kane et 
al. Hamiltonian, where instead the hole propagates by emitting or absorbing a single spin 
fluctuation. In our approach the spin is carried only by the boson a\ (which changes the 
spin by one), so that the conservation of Qz necessarily implies a quadratic Hamiltonian in 



the boson operators. A further simpliflcation of can be obtained following Ref. |15,16 



for a different but similar problem. One can exactly trace out the single fermion // from the 
Hamiltonian H using translation invariance and thus obtaining an effective spin Hamiltonian 
defined by the translation operator T^-^aiT^T-^ = ai^^-^ of spin-waves for nearest neighbour 
displacements 

^eff = ^ E(^ir.^o,., - mm.e^P^^ -J) + Hsw, (4) 

where i = denotes the origin and Hsw is the Heisenberg Hamiltonian that, at first order 
in 1/s reads: Hsw = £ J2 {'4'lj'^i,j - !)• 

Contrary to the undoped case, the Hamiltonian Hes cannot be solved analytically unless 
for the case t = 0, where H^s becomes quadratic [|l^,|l^]. However, a very good variational 
wavefunction that is exact in this limit, and which preserves all the symmetries of the 
Hamiltonian, is very easy to write down, in the form of the most general ground-state of a 
quadratic Bogoliubov Hamiltonian: 

I'^h) = exp{WB^,,aU]}\N). (5) 

Bij is non zero only if i and j belong to different sublattices (to fulfil Qz = 0) and its Fourier 
transform does not contain the modes aX k = (0,0) and k = (vr, tt) (to fulfil Q = 0) In 
order to evaluate and then minimize the expectation value of the Hamiltonian H^s over the 
state (H), one needs to evaluate both the average {'^ hl'iplj'ipijl'^ h) , and the average of the 
quadratic form ipljipij over |\E'/i) and the state |\E'^) = Tj-^|\E'/j) generated by the translation 
operator entering the kinetic term, i.e., {'^h\i'iji'i,j\'^h)- The state |\E'^) is clearly obtained 
by replacing in Eq. (^ the matrix Bij with B-'j = and because both the states 

l^'^) and l^^h) are gaussian, it is then possible to evaluate the averages in closed form. In fact, 
given two gaussian states {"^a) and l^&c) of the form (^), with A and C the corresponding 
matrices entering the exponential, we have 













c) 
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where 



(^^1^^) = det"i/2(/ - CA*) (6) 

Using the above equations is it easy to work out an exphcit expression for the expectation 
value of the Hamiltonian by simple linear algebra operations over the symmetric complex 
matrix Bi j. We have then obtained the optimal matrix Bij by minimizing the energy of 
the effective Hamiltonian ifgfj with the standard conjugate gradient technique. 

In table I we show data for the one hole energy (referenced to the undoped energy) 
and the quasiparticle weight Z = where I^^h) is the ground-state of Hsw- 

Because I^^h) is easily written in the gaussian form (|D [§,0, Eq. (||) allows to evaluate Z 
straightforwardly. The agreement of the spin-wave estimates with the exact diagonalization 
results [0 is surprisingly good, yielding a robust evidence of a finite value for Z in the 
static t = limit. 

The accuracy of the method remains very good even for t > 0, as it is shown in Fig.l 
for the quasiparticle weight at p = (0,0) and p = (7r,7r). We also see in Fig.l (b) a clear 
transition of the quasiparticle weight for the (vr, vr) momentum. Its value changes of about 
two order of magnitudes also in the exact diagonalization data. 

Our spin-wave approximation agrees with the exact diagonalization even in the details 
for t/J < 1. In fact, at the value t/J ~ 0.5, where in our simulation we find a singular 
point [see Fig.l (b)], there is a true level crossing in the exact diagonalization. The true 
ground-state is actually orthogonal (with different symmetry) to {"^h)- For larger t/J our 
approximate solution predicts that ^(^^7^) vanishes at a critical point for any finite size M. 
At the moment a similar analysis [|r^ for momenta close to (f , f ) (which is found to be 
always the ground-state for M —>■ 00, consistent with the general believe in the physical 
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region J not too small) indicates that Z^^il^il) remains finite up to t/J ~ 3 |T5|. However 
we cannot exclude a transition to d = at larger values oi t/J where our variational 
approach becomes unreliable. 

The fact that some of the low energy excitations of the Mott insulator may have a "non 
trivial" character may lead to a completely new classification of the charge excitations in the 
low doping regime. This surely requires further analysis and more analytical work. For the 
time being we point that a drastic change of the weight for momenta differing by the nesting 
wavevector g,r = (7i",vr), which are degenerate in energy, is a remarkable prediction of our 
approach which explains very well the numerical data on small systems and can be easily 
detected experimentally by photoemission experiments. For instance it is not possible to 
obtain the above property within the Kane et al. approach, because in this case the Green's 
function satisfies to G{k + q-j^iUj) = G{k,uj), so that the weights Zj. and Z^^g^ can only be 
equal. 

Work is in progress for extending the calculation at finite density of holes and/or at 
smaller values of J/t. 
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TABLES 

TABLE I. Linear spin-wave estimates of the quasiparticle weights and the one-hole energies in 

the static limit {t = 0) for the various clusters. The percentage relative errors refer to the spin-wave 
data as compared to the exact diagonalization results. 



M 


QP Weight 


% Error 


Enerffv 


%Error 


8 


1 


0.0 


2.5 


0.0 


10 


0.9919 


0.0047 


2.4308 


-0.35 


16 


0.9724 


-0.22 


2.3271 


-0.61 


18 


0.9688 


0.055 


2.3052 


-0.60 


20 


0.9637 


0.050 


2.2908 


-0.65 


26 


0.9512 


0.077 


2.2609 




32 


0.9422 




2.2415 






0.820 




2.17 
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FIGURES 

FIG. 1. (a) Quasiparticle weight as a function of t/J. The empty dots correspond to our 
spin-wave results; the full dots are the exact results on finite lattices. Continuous lines arc guide 
to the eyes, (b) Same notation as in (a) except that the empty dots refer to the 100 sites lattice. 
The dotted line is a local unstable minima for i/J > 0.6 
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